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Abstract 

The equilibrium susceptibility of uniaxial paramagnets is studied in a 
unified framework which permits to connect traditional results of the the- 
ory of quantum paramagnets, S = 1/2, 1, 3/2, . . . , with molecular mag- 
netic clusters, S ~ 5, 10, 20, all the way up, S = 30, 50, 100, ... to the 
theory of classical superparamagnets. This is done using standard tools 
of quantum statistical mechanics and linear response theory (the Kubo 
correlator formalism). Several features of the temperature dependence of 
the susceptibility curves (crossovers, peaks, deviations from Curie law) are 
studied and their scalings with S identified and characterized. Both the 
longitudinal and transverse susceptibilities are discussed, as well as the 
response of the ensemble with anisotropy axes oriented at random. For 
the latter case a simple approximate formula is derived too, and its range 
of validity assessed, so it could be used in modelization of experiments. 

1 introduction 

Modern magnetism, quantum mechanics and statistical mechanics have a long 
history of common development since the beginning of the XX century. Behind 
it are the names of the pioneers: Langevin, Brillouin, Bohr & Van Leeuwen, 
Landau, Van Vleck, Pauli, . . . After spectroscopy, magnetism helped the most 
in the understanding of the atom. Besides, some of the earliest successes of the 
Gibbs-Boltzmann distribution were indeed applications to magnetic problems 
[T]. Initially, applications to the simplest magnetic systems, atoms or small 
molecules with a permanent magnetic moment: paramagnets. 

When the magnetic moments are brought together they can interact (by 
dipole-dipole coupling or by "exchange" , due to the Pauli principle). Then a new 
game starts, with the possibihty of long-range order, spin waves, etc. Another 
important player, when the spins are placed in a molecular complex or in a solid, 
is the magnetic anisotropy. The spin-orbit coupling allows a magnetic moment 
to sense the electric field from neighboring ions. This gives rise to preferred 
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orientations, according to the space symmetry of the compound/soHd: cubic, 
tetragonal (biaxial), hexagonal (^uniaxial), etc [2 [3]. 

The magnetic anisotropy can also interplay with the spin-spin interactions, 
modifying the spin-wave dispersion relations or determining the spin orientation 
inside ordered "domains" |H [3]. If there is no room for domain walls, as in 
sufficiently fine particles, the anisotropy would dictate the stable orientations 
of the whole magnetic moment (the dipolc-dipole interaction also plays a role, 
as it "adds" to the anisotropy). Those stable orientations are the basis of using 
such particles as physical bits. 

1.1 quantum &: classical superparamagnets 

One more element sets in if the temperature is high enough (but still below the 
ordering temperatures), or when the particle is sufficiently small (~nm). Then, 
by thermal activation, the net moment can overcome the anisotropy energy 
barriers and flip back and forth between the minima [H [5] . This is a nightmare 
for magnetic recording [6], but a blessing for new phenomenology. Inside we 
have an ordered magnet, which seen from the outside resembles a paramagnet, 
but with a very large magnetic moment, S ~ 100-1000 — superparamagnetism. 
And yet another connection with statistical mechanics is established, as the non- 
equilibrium orientational distribution of those "classical" spins is described with 
the tools of Brownian/stochastic dynamics, e.g., Fokker-Planck or Langcvin 
equations [3111111. 

The 80s and 90s brought, through the advances in chemical synthesis, a new 
member of the paramagnetic family, molecular magnetic clusters [10| (also rcf- 
ereed too as single-molecule magnets). These are made of complex molecules 
with a net spin 5* ~ 5, 10, 20, somewhere in between the traditional para- 
magnets and the superparamagnets. But, contrary to magnetic nanoparticles, 
they are assembled in molecular crystals, minimizing extrinsic sources of disper- 
sion (in size, anisotropy parameters, etc.) and turning the comparison between 
experiment and theory cleaner. 

As for the static part of the magnetic anisotropy, it can be quite large in 
molecular magnets, so the stable orientations could provide the basis of a quan- 
tum bit. Nevertheless, environmental effects (decohercncc in particular), emerge 
already from their own nuclear spins or from the dynamical part of the spin-orbit 
coupling (neighboring ions are never static; the lattice oscillations modulate the 
local electric field, which modifies the orbital motions, eventually affecting S). 
Again, a curse for technology, but a blessing for physics [11]. The total system 
can be described by a spin-phonon Hamiltonian [12l [131 [El j a-nd constitutes a 
neat realization of the spin-boson paradigm of the theory of quantum dissipa- 
tive systems |15| : but with a S > 1/2 immersed in a bath of true "harmonic 
oscillators" (the lattice phonons). 
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1.2 the spin Hamiltonian 

The theoretical description of paramagnets starts with a spin Hamiltonian, 
where one replaces the true microscopic Hamiltonian (or a part of it) by one 
written in terms of spin energy levels and operators (a procedure initiated by 
Heisenberg, Dirac & Van Vleck [ini § 72]). This provides an equivalent descrip- 
tion of the statics (based on the spectrum) and of part of the dynamics. The 
inner complexities, symmctrization of wave functions, local fields, spin-orbit 
coupling, are left in the backstage and we only sec their effect through the spin 
Hamiltonian [T1IT2l[T6]. 

For many classical superparamagncts [l7j and molecular clusters [TU] , a good 
first approximation is provided by 

■Ho = -DSl-B-S , (1.1) 

with D the uniaxial anisotropy constant. D > Q produces two preferred spin 
orientations (the "bit" , Fig.[T]) and brings in the physics of bistable systems (e.g., 
thermal activation over a potential barrier). But there is also a rich "coherent" 
or transverse dynamics [TH1|TH], with spin resonances, absorption peaks, etc0 

1.3 the magnetic susceptibility S>c linear-response theory 

Most of what we know about magnets, and even their names and the ways we 
classify them, is based on their characteristic response |T]. That is, on how they 
respond to the application of controlled magnetic field probes. When the probes 
are small enough, linearizations and perturbation theory, in its different guises, 
can be used by the theorist. 

In molecular magnetism a procedure due to Van Vleck [U Ch. 3] became 
quite popular. The energy levels of a paramagnet, in the presence of the prob- 
ing field, are obtained using quantum-mechanical perturbation theory, and the 
result plugged into a Gibbs-Boltzmann formula. In other, more complicated 
problems of theoretical magnetism (e.g., quantum spin chains [HI [231 HH US] or 
inner structure of molecular clusters |26| ) other tools from linear response the- 
ory are employed. The Kubo correlator, in particular, provides a formal general 
solution for the response functions, the susceptibility x, without invoking per- 
turbed levels at all (the usual "magic" of linear- response theory [?7l I28j ). 

1.4 our approach and goals 

We aim to calculate the equilibrium linear response for paramagnets with simple 
anisotropics. But in a way that could connect the traditional results for quantum 
paramagnets, all the way up to the results of classical superparamagnetic theory. 

^ The big spin Hamiltonian 1 11. Il l is also an effective description of a collection of N two-level 
systems (not necessarily magnetic), used in atom optics since the 70s 1191 . Non-interacting 
entities correspond to D = 0, while D with a transverse field Bx describes certain types 
of uniform interaction (the Lipkin-Meshkov-Glick model originally from nuclear physics). The 
theoretical work on this front was revived with the boom of cold gases and condensates I20II21| . 
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Figure 1: Energy levels £,„ ~ —Dm? — Bzm of anisotropic spins S ~ 1, 3/2, 2, 
and 4 at zero field. For 5=1 and 3/2 we set D = 1. For S = 2 and S = A, the 
anisotropy constant is scaled as D = so that these spins have the same 

energy barrier as = 1. 

The approach of Van Vlcck, at least in its popular form of diagonalizc & lin- 
earize [29l [30], becomes awkward for moderately small S. Therefore wc will 
resort to the Kubo correlator formalism, and particularize it with the (su- 
per)paramagnctic problem in mind. 

Schematically, the approach takes two steps (Sec. [5]): (1) Obtain the cor- 
relator Jgda Tr [e*^^^"")^" V e+°'^° V] , and sandwich here and there identities 
I = '^m\m){m\ formed with the unperturbed levels Ti.o\m) = e,„|m). This 
gives a susceptibility x ~ n c"'^^'" |KinP^[/?(£m - £n)] with the Kubo 

function K{X) = {e^ — 1)/X evaluated at the level differences. (2) Next, for 
uniaxial problems, Ti.Q — ~D Si ~ Sz, the matrix elements of the perturba- 
tion V simplify, Vmn ^ y\\^m,n + V^5ra,n±i, and the above x can be reduced to 
a single sum over the unperturbed spectrum (in the line of naive linear-response 
theory). One component of Xi the longitudinal susceptibility, comes from the 
part of V commuting with Hq; the non-commuting part gives the transverse 
susceptibility [Eqs. ^T^i^U^]- 

Some advantages of this approach are: (1) The intermediate step of getting 
eigenvalues in presence of the perturbation is bypassed. (2) There is no small- 
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denominators problem; K{0) = 1 ensures automatically the finiteness of the 
perturbative result, without the need of manual pairing of degenerate levels or 
the like [SI]. (3) The susceptibility formula, expressed as a single sum over the 
spin levels m, is ready to use for any S. 

As a check, we recover in Sec. |3] a number of traditional results for param- 
agnets, including the transverse response of small spins S* = 1 and S ~ 3/2. 
Then, we will move on to arbitrary 5* in Sec. |4] and discuss some features of 
the temperature-dependent susceptibility curves, like crossovers, peaks, and the 
isotropy of x in the randomly oriented ensemble. We will study how these 
features evolve with S under appropriate scalings (Fig.[T|). 

The interest in the x(^) curves stems from their routine use as indicators 
of anisotropies, interactions between spins, etc. [29l [30l |32] . Therefore it is 
important to understand well the paramagnetic "ideal gas" (as the implicit 
reference curve in those procedures), including the effect of simple anisotropies 
(e.g., uniaxial, biaxial, etc.). We hope that our results will convince the reader 
of the advantage of having a unified framework for the response of the whole 
paramagnetic family, from the traditional members, then molecular magnetic 
clusters, up to classical superparamagncts. 

2 derivation of the susceptibility 

Instead of starting writing down directly the linear susceptibility in terms of a 
Kubo correlator, we will briefly discuss its basis on plain statistical-mechanical 
perturbation theory. This will make the presentation more self-contained, and 
will hint at how to proceed if higher orders are required, beyond linear response. 

2.1 statistical mechanical preliminaries 

Statistical mechanics is about computing averages, variances/dispersions, etc. 
for a given system with the sole input of the system Hamiltonian TL and a 
few parameters incorporating the effect /presence of a thcrmalizing environment 
(like /3 = 1/T, chemical potentials, etc.). The averages are generated from the 
density matrix g oc e~^'^, and for a given operator quantity Q one forms the 
following basis-independent object 

{Q) =Tr{gQ) ^Y^QmnQurn ■ (2.1) 
ran 

For consistency of the definition, g is "normalized" to have unit trace Tr (g) = 1 
(so the average of a constant is the constant itself): 

g = e-^'^/Z, Z = Txc-^^. (2.2) 

The normalization function, Z, plays a central role in several parts of the for- 
malism. In particular, as an auxiliary/ abridged generator of averages (by dif- 
ferentiation) with the advantage of being a scalar instead of an operator. 
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2.1.1 response observables 

Let us imagine that the quantity we are interested in can be obtained from the 
Hamiltonian by differentiation 

Q = —dfH (operator identity) . (2-3) 

Say, the total Hamiltonian is of the form Ti. = Ti.Q — / Q, as in the examples of 
force-coordinate, —Fx, or ficld-dipolc/spin, —B-S. Then the average of Q can 
be obtained from the partition function as follows 

(Q) = |^Tr(a,e-n = |a(,/)Tr(e-^«) = | . (2.4) 

Here ( )' denotes derivative with respect to the "thermal" force/field parameter 
^ := 13 f entering in = -PHo + ^Q.E 

We will refer to the average of a quantity Q derivable from the Hamiltonian 
as above as a response observable. This definition does not include pure thermal 
quantities, like the thermodynamical energy or the specific heat, although these 
can be obtained by differentiating Z with respect to /3 = 1/T. 



2.1.2 susceptibility & derivatives of Z 

As the force/field parameter / is at our disposal, we can change / and from the 
induced change in (Q) learn how the system "responds" . A natural quantity 
then arises 

X = Of iQ) , (2.5) 

quantifying the sensitivity of the system to changes in /, and for this reason 
called the susceptibility, x is thus a second derivative of Z, and if we keep on 
using £, = (if and ()' = 9^, we have x = I ^) ~ I ^Y] - where the square 

brackets is merely [Z' / Z)' . 

We can think that / is the part of the force/field that we change (the probe), 
over a fixed bias /o. Then, from — (/o + f)Q we can move —foQ into the 
unperturbed part Ho, and eventually evaluate our derivatives at / = (linear 
susceptibility). Then, all we need to do is to obtain the ^-expansion of Z to 
second order, Z ~ Zq + ^Z\ + ^^Zi, and the linear susceptibility will simply 
follow as a combination of the expansion coefficients 



~Z^ 



£2 



£1 



(2.6) 



^ The proof of (12.41 1 is a bit less direct than it looks like, because one cannot differentiate 
the exponential of a matrix plainly as 9j c^'^' ^ C^'^'S^yi, when A and d^^A do not com- 
mute. However, using the proper identity 9^ c^'^' = f^Aa (9^ A)e+'^'^(^\ with a 
an auxiliary parameter, along with the cyclic property of the trace Tr (C -D) = Tr (D C), one 
arrives at the desired result; (9^Tre'^'^' = Tr (e^9^j4), with the trace removing any ordering 
problem. 



6 



2.2 perturbative treatment 

When one proceeds to expand the partition function 

Z = Tr c-^'(^°+^) , V = -fQ , (2.7) 

with respect to ^ = f3f, one faces the problem of handling exponentials of 
operators or matrices. Classically, one proceeds by factorizing and expanding 
ga+6 = e" e& 2± e'^(l + 6 + + • • • ). For operators, however, plain factorization 
does not hold, and one has instead a Baker-Campbell-Hausdorff formula 

^A+B^^A^B^-:^[A,B]^+^aAAA.B]] + [BAA.B]])_^ . (2.8) 

This kind of expression is useful if one can recurrently write higher order commu- 
tators, e.g. [A, [A,B]], in terms of lower order ones. A, B, [A,B]. For example, 
when A ^ Sz and B ^ S± = Sx i iSy (isotropic spin), or when A ^ b+b^ 
and B ~ 6±, with [fe„,6+] = I (harmonic oscillators). But in the general case 
equation (|2.8p is of little use to do perturbation theory in one of the operators 
(actually, already for A ^ Si and B ^ S±). 

2.2.1 perturbations from a Kubo identity (interaction picture) 

A way out in the problem of handling operator exponentials is provided by Kubo 
identities of the type [171 p. 148] 



e^+^ = e^ 



dae-'^^Be+"(^+^) 







(2.9) 



where one resorts to an integral over an auxiliary parameter|f| The integrand 
includes B (now freed) and again the exponential of the sum. Therefore, by 
iterating the same expression, e<^iA+B) ^ ^aA _|_ J^'^dr c^^-^ B c+^i^+B)] one 
generates the successive powers of B and can do perturbation theory. 

For example, to get second-order derivatives of Z, one just needs to iterate 
to second order. Then the B in the last e+'^(^+^) is dropped (as there are 
already two B multiplying it), and one is left with 



^A+B ^ ^A 



I + / da c-"^ Bc^"^ 



• 1 



+ / dac-"^ Bc+"^ / dre-"^Be+"^ 



(2.10) 



Here we recognize the "interaction-picture" evolution e~'^^ B e^"^^ . Indeed the 
structure of the above formula is present in most of quantum mechanical pertur- 
bation theory [33], including scattering; it is also used to derive weak-coupling 
master equations in quantum open systems [^ [33] . 

^ The same kind of auxiliary parameter, with in principle no physical meaning, we used 
to handle 9je'*(^'. The proof of Eq. 112.91 1 follows in four simple steps: (i) Isolate the in- 
tegral by left multiplying by c'^. (ii) Define F{\) := e"^ e+^(^+S) ; then F{0) = I, 
while -F(l) is the target, (iii) Differentiate with respect to A (no ordering problems) to get 
wc arc almost done [cf. the integrand in Eq. ||2.9|| ]. (iv) Integral 
reconstruction -F(l) = F(0) + /o^dcr F'(u) gives the [l + ■ ■ ■ ] on the right-hand side of l|2.9|l . 
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2.2.2 tracing & the Kubo correlator 

Tracing is now required to get Z = Trc"''^. The trace will simplify the per- 
turbative treatment (as compared with perturbed time evolutions), by allowing 
to move operators around in combinations like e~'^^ B e'^'^^ , using the trace's 
cyclic property. We will undertake this first, arriving at the structure of the 
Kubo correlator |35| , and then simplify further by doing the trace in the eigen- 
basis of A. 

i) simplifying the trace by cycling. Using Tr (C D) = Tr {D C) when 
tracing the second term in (|2.10[) one gets rid of the integral, as the first ex- 
ponential can be moved to the end, canceling the dependence on the auxiliary 
variable a. The third term also simplifies, following a procedure explained in 
most field-theory books Q Collecting the results one arrives at 

Tr(e^+^) ^ Tr (c-^) + Tr (c^ B) + i /'da Tr [e^^-")^ B e+"^ B ] . (2.11) 

We see that one of the auxiliary parameters, a, is still with us; indeed the last 
term is a bare form of the Kubo correlator, /g'^dcr Tr [c'^'^"'^)^" V c+'^^o 1/] . 

The goal is therefore accomplished: tracing the exponential of a sum of 
operators, to second order in one of them. For this Bakcr-Campbell-Hausdorff is 
not suited, while a classical handling of the operators 0-^+^ 2± c^(1 + B + BV2), 
would give correct results only to first order. 

ii) tracing in the unperturbed eigenbasis. When the eigenstructure 
of the operator A is known, A\m) = am]™), the trace can be written explicitly 
in terms of the eigenvalues and the matrix elements of the perturbation 

A\m) ^ ara\m) Bmn ^ {m\B\n) . (2-12) 

Indeed sandwiching identities I ~ between e^ and B, the three parts 

of Eq. (pHI) give 

Tr (e^) = 51 e"^" > Tr (e^^ = (^.13) 

7n rn 

\ /''daTr[c(i-^)^Bc+'^^B] = i ^ e"" |i3„„|2^(a„ - a„) , 

where we have written B^nnB^m ~ l^^mnl^ for Hcrmitian B. The integral over 
a produced the Kubo function 

K{X):= [ dac'^^ ^ {c^ -1)/X , (2.14) 

* See for example |36l Ch. 8.3]. One first converts the integral over the triangle a € [0, 1] 
T € [0,0-] into the integral over the unit square [0, 1] X [0, 1], by using "chronological" ordering 
(this yields a factor 1/2). Then the cyclic property of the trace is used to show that the 
integrand does not depend on one of the integration parameters. 
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which enters evaluated at all eigenvalue differences X = a„ — a„i. This function 
will follow us all the way to the final expressionsH 

2.2.3 application to the original perturbative problem 

Let us write down the explicit correspondence with our original perturbative 
trace problem Z = Ti exp[— /3(7io ~ / Q)] 

- (3n^ ~(3no+^ Q . (2.15) 

A B 

Now the unperturbed eigenstructure reads 

Ho|m) = emlm) Um = -PSm Bjnn^CQmn- (2-16) 

Then comparison of Eqs. (|2.1ip and (|2.13p with the ^-expansion of Z ~ Zo + 
^Zi + i^^Z2, gives the sought for coefficients 

m m m n 

(2.17) 

In the last term we have introduced the shorthand 

Kmn ■■= K[P{em ~ £„)] A-(/3A„„) , (2.18) 

with the Kubo function evaluated at the level differences A^n '■= £m ~ £n 
("transition" frequencies). 

i) absence of small denominators problems. A final remark on per- 
turbation theory with close levels. If we combine the denominator — e„ 
from Kmn with the matrix element IQmnP; one actually sees the structure of 
plain quantum-mechanical perturbation theory: QmnQnm/ [sm — £«)■ How- 
ever, there is no need to handle degenerate levels, if existing, in a special way 
(i.e., no need of pairing, etc. [31]). The formalism ensures that the finite tem- 
perature perturbative treatment is finite as well. Indeed, Taylor expansion of 
K{X) = (e^ - \)IX 

K{X) ~ 1 + X/2 + X'^/6 + • • • , (2.19) 

shows that degenerate levels A,„„ = would contribute a finite A'(0) = 1. This 
property spares us with a degenerate perturbation theory to handle close levels 
Sm — En, and is built in the finite-temperature formalism (to second order at 
least). 



5 Note the connection K{X) = l/Wi(X), with Wk(X) = X*/(e^ - 1) a transition rate 
from the theory of open quantum systems (fc = 1 corresponds to an "Ohmic" bath 1151 1341 '!. 
Thus, the Kubo function enjoys a "detailed balance" relation as well: K(—X) = e~^ K{X), 
useful in some manipulations. 
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2.3 formulas for the susceptibility 

Now it is immediate to write down explicit expressions for the susceptibility, by 
plugging the Z expansion coefficients (|2.17p into x = /3[(-22/-Zo) — (-^^iZ-Z^o)^] • 
This form is generic for what we called "response observables" in Sec. 12.1.11 
that is X = <9/ (Q) with Q = —dfH, as in the examples of a coordinate, —Fx, 
a spin component, —B^Sz, etc. 

In what follows wc will address the structure of the resulting susceptibilities, 
and then particularize the discussion to "ladder perturbations" , which include 
simple mechanical oscillators and uniaxial paramagnets. 

2.3.1 general x zero bias 

As a quick illustration, in the unbiased case where F = or i?^ = 0, one has 
2^1 = as well, and the susceptibility simply reads (we will restore Zi ^ 
shortly) 

X= ^ Ve-^-'" |g™„pi^„„ . (2.20) 

m n 

This is the kind of ready-to-use expression we mentioned in the introduction. It 
is written fully in terms of the unperturbed eigenstructure {Sm, Qmn = 

{m\Q\n), and K^n = K[(3{em — e„)]. It looks like a average over the unper- 
turbed system, resembling the classical result x ~ {'^i)^ but with Kmn encod- 
ing effects of non-commutativity of the perturbation Q and the base Hamil- 
tonian. With simple numerical diagonalization, the x above can be used for 
a non-linear oscillator, "Hq ^ —^kx^ -\- \qx^, or arbitrary anisotropic spins, 
Tio ~ ~D S^,^ — ± - ■ ■ , with S± = Sx ^ iSy, in particular biaxial systems. 

2.3.2 X foi" ladder perturbations 

(harmonic oscillators and uniaxial magnets) 

The general expression for the susceptibility obtained by plugging the Zi from 
Eq. (|2.17p into x = P[{22/Zq) — (Zi/Zq)^] , simplifies when the coupled observ- 
able Q gives, at most, transitions to adjacent levels: m ^ m and to — > to ± 1. 
Then one can write 

Q = b-S^boSo + ^{b+S-+b-S+) , (2.21) 

with the following action on the unperturbed basis 

So\m)^e°Jm) S±\m) =e^\m±l) . (2.22) 

The harmonic oscillator corresponds to no central term £^ = and the creation- 
destruction factors £^ = [(to -|- i) ± ^]^^^. For spin problems, the coeffi- 
cients are i'^ = m and the custom angular-momentum ladder factors = 
^S{S+l)-m{m±l).E 

^ Note that here wc arc already thinking of an unperturbed Tio = Ti-oiSz), otherwise the 
standard basis, where S has the ladder properties 1 12. 221 1. does not diagonalize Tio, and we 
would have to resort back to the more general 1 12. 20t . 
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For these "ladder perturbations" the terms with bo in Z2/Z0 and (Z2I Zq)^ 
can be combined, while the ladder action (|2.22p reduces the double sum X^mn 
to a single sum. The susceptibility then reads 

X= bl A^e-''-[(Cf-(C)^] 
+ i6+6_ ^Y.{[ll;)^G^^^^^ + [l-;)'Grn^^.{\ (2.23) 

^0 

where we have introduced the shorthand Gm^n = c~^^"^ K„in- It is not dif- 
ficult to check the symmetry Gn,m — Gm,m from the detailed-balance prop- 
erty K{~X) = e~"^ K(X) mentioned above. This symmetry, together with 
— ^m±i: leads to the two sums in the transverse 6+6_ part being equal (the 
sums, not the summands). Therefore we can keep one of them, replacing the 
factor 1/4 in front by 1/2.0 



2.3.3 final expression/summary for uniaxial magnets 

We conclude writing explicitly x for paramagnets with TCq = HoiSz) probed by 
a field oc ~b-S. The unperturbed basis is then the standard basis Sz\m) = m\m). 
Besides 6|| := bz is the direction cosine of the probing field parallel to the magnet 
local axis, while b\ ~ b^b^ corresponds to the transverse one. 
The susceptibility can then be decomposed into 

X^b\x\\+blx^., (2.24) 
with the longitudinal & transverse components given by 

XII 

Here Zq — ^-0^^ gj^^ have done the mentioned reduction of the dupli- 
cated terms (keeping X]m(^m)^^m,m+i) ^"^^ simplified some notations 

4 := (^+)' = ^(^+l)-m(m + l) (2.27) 
■■= /v,„,,n+i =X(/3A™„+i) . (2.28) 

Thus the Kubo function K{X) = (e^ — 1)/X enters evaluated at the "upward" 
transition frequency l^mm+i = Em ~ £m+i- The rest of the article will consist 
essentially of examples and application of Eqs. (|2.25p and (|2.26p . 

^ The equality J]m(^m)^Gm,m+i = Z],„(^m)^Gm,m-i can also be proved directly from 
the Kubo correlator, as those terms correspond to j'^dalr [c(P-'^)'^o(S^) (,+crno(S^) g_^y 



P. 

JL 

2Zo 



-Per, 



F K 



(2.25) 
(2.26) 
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some remarks. Before closing this section, two remarks are in order. First, 
the longitudinal x\\ can be obtained directly from derivatives of the unperturbed 
partition function Zq, with respect to the static field, as it is well known. In 
our notation, when TLq = HaiSz) — Sz B^, the levels are —fiem = —P^m + my, 
with the dimensionless variable y = [3Bz- Then Zq = X^m whence 
Zq (m^) = d^Zo/Ay^ for any moment (m^). As a result 

, , 1 AZn I o> , ,9 d (to) 
{m)^—-^ {^^^^) - {m)- = ^ , 2.29 
Zq dy dy 

from which x\\ follows as X|| = f3d (m) /dy with y = PBz. 

The last remark is on angular behavior. The susceptibility proper is a tensor 
quantity Xij relating two vectors (magnetic moment and probing field). What 
we have been using throughout is the projected scalar form x Xij^ibj = 

Xzzb^ + Xxxib^ + by), where the last form holds for uniaxial symmetry. Then, 
for a system of non-interacting spins with a distribution of axes orientations, 
the scalar x follows from the corresponding angular averages 6| and b\. This 
can be done, for example, when there is no bias field singling out a preferred 
direction; then one has 

X = iX|| + |x±, (2.30) 

for the susceptibility of a system with anisotropy axes distributed at random 
(powder sample, liquids, etc.). 



3 paramagnets 

In this section we will check the formulas of Sec. [5] particularizing them to 
simple paramagnets. We will consider isotropic spins VS" and anisotropic spins 
with = 1 and S ^ 3/2. The susceptibility for anisotropic problems is typically 
obtained using Van Vlcck's method, solving the eigenvalue problem in presence 
of the probing field i30j . We will plainly recover those standard formulas 
of molecular magnetism textbooks (bypassing diagonalization) and minimally 
extend them by including longitudinal bias fields. Here we will also see the 
first examples of the x(r) phenomenology that will be discussed later on for 
general S. □ 



3.1 isotropic spin (Brillouin x) 

This is the simplest paramagnetic problem [1] [T^], with Hamiltonian TYq = 
—Sz Bz and spectrum —fiem = my. The energy levels are equispaced /9A„i„i+i = 

* The reader may well wonder why we have not used Bosc/Fcrmi statistics for integer/half- 
integer S. First, the only basic statistics is Gibbs— Boltzmann, with Bose/Fermi distributions 
as particular/worked cases of oc e~l^'^ for the ideal gas of indistinguishable bosons or 
fermions. But we implicitly bypass indistinguishability by locating/labeling each spin at a 
given lattice point (localized moments' magnetism) [T] Ch. 3.1] |12l Ch. 1.2]. This is the same 
approximation routinely used in solid state physics and quantum chemistry, where one does not 
symmctrizc/antisymmctrize the states with respect to the nuclei exchange; an approximation 
grounded on the non-overlapping of the ions' wave-functions for sufficiently localized states. 
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y, yielding a m-independent Kubo factor K{y). Besides, the partition function 

^0= E e- = ^^%^#^^ y = B^iT, (3.1) 

m— — S ^ A ^ ' 

follows readily by summing the geometric series X]m(^^)™- 

3.1.1 longitudinal susceptibility 

As discussed above, the longitudinal response follows entirely from -derivatives 
of Zq. The first moment reads (magnetization) 

(m) = (^ + i) cth \{S +\)y\-\ cth (i y) , (3.2) 

with the right-hand side defining the Brillouin function. One further derivative, 
using (cthx)' = 1 — cth^x, gives the susceptibility X|| = /9d (m) /dy as 

X||-/3^(^+l)-/3[(5+i)'cth2[(5+i)jy]-icth2(iy)] . (3.3) 

In the weak-field limit, one can use cth^x ~ 1/x^ + 2/3 to show that X|| duly 
reduces to the Curie law 

XII Xc = ^l3S{S+l) . (3.4) 

This famous 1/T dependence was found experimentally by Pierre Curie in 1895 
[1] Ch. 3] and it has been fruitfully exploited for calibration and thermometry 
in the low temperature world. It expresses the decrease of the response with 
increasing T due to the thermal misalignment of the dipole moments away from 
the probing field direction. 

3.1.2 transverse susceptibility 

From the general Eq. (|2.26p for x±, plus the equispaced —fiem = my for 
isotropic spins, one forms 

x^ = f E— (3-5) 

The m-independent Kubo term K{y) ~ (e^ — 1)/^ can be taken out of the sum, 
which can be done exphcitly: 

y 

E -r- [^(^ + 1) ~ "^(™ + 1)] = [ct^ (5 y) - 1] (™) • (3-6) 

rn ^" 

But the square bracket can be written as 2/(e^ — 1), which combined with the 
Kubo factor leaves the simple form xi- = {P/v) ('^)- That is (recall y = (33^) 

X^ = ^ '-^ lPSiS + 1), (3.7) 
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which also recovers the correct Curie Umit at low fields, as expected from the 
restoration of isotropy@ 

It is important to remark that x± docs not follow from the transverse fluc- 
tuations of the spin, as one might naively expect. Indeed, using S*^ + Sy = 
S{S + 1) - S*^ one gets for Ho = -S^ 

13{SI) = 13{SI) = \ V cth (1 2;) X ^ ^ XX . (3.8) 

They only agree at high temperature, y ~ (iBz «C 1, where the extra prefactor 
i y cth (i y) — > 1. However, at low temperature, x± ^ S/B^ (following satura- 
tion of the magnetization (m) — > S), whereas (3{S'^) — + S/2T, which can grow 
without bound as T — *■ 0. Two different behaviors indeedF°l 

3.2 anisotropic spins 

When a paramagnetic ion is embedded in a molecule or a solid, the spin finds 
preferred orientations which depend on the symmetries of its neighborhood. 
This magnetic anisotropy can be described by adding to the spin Hamiltonian 
a term Hs.{S) with "reflection" symmetry Hs.{~S) = Ha (5) (to respect time- 
reversal invariance). The simplest model is the uniaxial Hamiltonian 

no = 'DS^,-S,B,, (3.9) 

with D the anisotropy constant. To do statistical mechanics we introduce the 
dimless d = D/T k, y = Bz/T, and write —Psm = dm? + ym. 

The minimal mathematical extension from linear to quadratic in S has im- 
portant consequences. First, the spectrum Em = —Dm? — B^m can have a 
single well (D < 0, "easy-plane" anisotropy), or it can display a bistable struc- 
ture {D > 0, "easy-axis"; see Fig. [T|). Second, the energy levels are no longer 
equispaced, Amm+i = D(2m, + 1) + B^., becoming closer near m = (no "har- 
monic oscillator" equispaced simplicity). In dynamics this gives a multiplicity of 
precession frequencies (absorption peaks), tunnel splittings and relaxation rates 
[T2I [m [37l [38] . But already in the statics we will find a m-dependent Kubo 
factor i^(/3Amm+i): expressing that e~''^' S±c'^'^^~ cannot be written in terms 
of S± only (another reason underlying the simplicity of isotropic spins). Finally, 
the partition function Zq = EmC'''"'^^™, with d = D/T and y = B^/T, can- 
not be summed explicitly (Gauss type sums), though simple formulas can be 
produced for small S. 

^ For S = 1/2 (the two-level system), one has (m) = i th(^y) (population difference), 
so that x\\ = 3/^/ ch^(^ J/), while x± = th(^y)/2Sz. Both duly give Xc = at zero 
field [S(S + 1) = 1/4], In the opposite, classical limit S ^ 1, one has X|| = PiJ?L'{£) 
and XX = [t^/ fo'" isotropic superparamagnets, with the Langevin magnetization 
L(g) = cth^ - 1/^ and g = fiB^/T; see for instance [g] Eq. (3.74)]. 

The unbounded /3(5j) can be seen as due to "zero point fluctuations", since + 5^ — > 
S{S + 1) — S'^ = S, is different from zero even for "fully aligned" Sz —* S; the actual x± 
corrects for this, and leaves the response induced by the probing field. 

Slow dynamics at low T is an important consequence of the anisotropy. From the point 
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Table 1: Energy levels — /3em = dm"^ +ym, ladder factors £^ = 2 — m(m + 1), 
and transition frequencies f3A„im+i = /3(£m — £m+i) = d{2m + l) + y, for 5 = 1. 
Note Sis + 1) = 2 and lS{S + 1) = 2/3. 

3.2.1 susceptibilities for 5* = 1 

Spin one is the first case giving new phenomenology; in the case S* = 1/2 the 
term —D Si produces a uniform level shift, so that x is the same as that of the 
isotropic spin discussed above. 

For 5 = 1 at zero field, there are two facing energy levels m = ±1 (minima 
for D > 0) with the level m = in between (a maximum or minimum, depending 
on the sign of D; top left panel in Fig. [T]). In Table [1] we have collected the 
energy levels, ladder factors, and transition frequencies required to calculate 
the equilibrium response. In terms oi d = D jT and y = B^/T ^ the partition 
function (12.21) reads 



Zo = 1 + 2e'' chy , (5* = 1) , (3.10) 
which enters in both X|| and 

i) longitudinal susceptibility X\\{S — 1)- Differentiating Zq one gets 
the magnetization (m) = Z^^dZQ/dy, which can be cast in the suggestive form 

ch y + i e " 

The m — Q level does not contribute to the numerator ^ m e~'^^™ , but occu- 
pies a "phase space" ^e"''. When for d ^ 1 the level m = gets thermally 
depopulated, one recovers (m) ~ thy, as in a two-level system. 

Differentiating (m) gives the fluctuations {iii^^ — (m) = d (m) /dy, which 
times j3 = \ /T gives the longitudinal susceptibility 

1 + ic^'^chy , , ^ 

^"-^(Th^TW' (^-D- (3-12) 

Again, d ;3> 1 leads to the two-level type susceptibility X|| — P/ ch^y. 



of view of the many relaxation rates 1371 1381 . one of the rates is well separated from the 
others (the analogue of the lowest non- vanishing eigenvalue in the framework of the Fokker— 
Planck equation). The faster rates correspond to intrawell modes, while the slow mode is 
associated with the overbarrier dynamics of the spin; due to thermal activation it is suppressed 
exponentially at low T, by exp(— ACZ/T). The phenomenological "blocking temperature" is 
that where the observational time window tm matches this slow dynamics; below it, the given 
technique does not "record" anymore equilibrium properties. 
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ii) transverse susceptibility x±{S = !)• Picking from Table [1] plug- 
ging in Eq. (|2.26[) . and playing with K{—X) = c^^ K{X), we can write the 
transverse response as X-L ~ (/3/-2o) — j/) + K(d + y)] . Alternatively, we 
can unfold the Kubo fmictions and use {d±y)/l3 = D±Bz obtaining 



1 \c'^-y - 1 c'^+y - 1 



As mentioned before, X-L remains finite when adjacent levels become degenerate 
(e.g., when — ±D). Note finally that when d — *■ 0, we have 

X. - - ^ , (3.14) 

chy+ ^ Bz 

duly recovering the transverse susceptibility p.7p for isotropic 5* = 1 spins. 

iii) temperature dependence of x\\ ^^nd xj. for S — 1. Wc will set 
Bz = Q to describe briefly the features of the T-dependent susceptibility. The 
above expressions, Eqs. and ((Hl^ . then reduce to (cf. Ref. [IHl Eq. (2.9)]) 

2 1 2 1 - C^/^ 

^Il = r 2 + e-^/^ ' ^^ = :P 2 + c-^/^ ' (^ = ^)- (^-^^^ 

These susceptibilities arc plotted in Fig. [H At high temperature both go over 
the Curie curve Xc = <5'('5' + \)/'iT. This can be seen from the formulas too, 
as \D/T\ < 1 yields X|| - 2/3r and X-L ^ 2/3T (now 5(5 + 1) = 2). As the 
temperature decreases, X-L shows no peaks and tends to the constant response 
Xi_ = ^/D (the analogue to the classical torque susceptibility X-L ~ i^e^ 
The longitudinal response, on the other hand, goes over the 2-state asymptote 
X|| ~ 1/T = S'^/T at low temperatures. Therefore, as S'^ > S{S + l)/3, the 
increase of X|| is faster than 1/T in the intermediate range. 

We mentioned the routine use of x{T) curves as characterization tool [32] . 
For instance deviations from 1/T laws can indicate spin-spin interactions (think 
of the mean-field x ~ ^/(T—O) when observed over a short T interval). However, 
one should be careful with not overlooking other possible sources of deviation, 
as those due to the anisotropy just discussed. 

3.2.2 susceptibilities for S = 3/2 

Wc move on to the next spin value. The anisotropy term —D Si gives two pairs 
of degenerate levels, without a central maximum or minimum (see Fig. [l] top 
right). The pair |7ti| = 3/2 is above or below the pair |to| = 1/2 depending on 
the sign of the anisotropy constant D. 

We have put up in Table [2] the relevant quantities to calculate the suscepti- 
bilities of this 4-level system. From there we first compose the partition function 
{d = PD,y = f3Bz) 

Zo = 2[ch(2//2) + e2'' ch(3y/2)] , (5 = 3/2), (3.16) 

which will be needed in both the longitudinal and in the transverse response. 
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Figure 2: Susceptibilities vs. temperature of spins 5* = 1 (left panels) and 
S = 3/2 (right) with easy-axis anisotropy D = 1 (see the spectra in Fig. [1]). 
Top panels: Longitudinal and transverse susceptibility vs. T. The Curie law 
S{S + l)/3r (thin solid lines) is approached by both x\\ and x± at high T. At 
low temperature X|| goes over the 2-state behavior S'^/T (dashed thin lines). 
The transverse susceptibility tends to the constant x± — l/[2-D(l — 1/25)] at 
low T [Eq. (|4.8p ]. Bottom panel: Susceptibilities plotted as Tx, including the 
average for axes distributed at random, x = ^X\\ + f X-Lj which shows deviation 
from isotropy at lower temperatures. 



i) longitudinal X|[(S' = 3/2). The y-derivative of the above Zq produces 
the magnetization 

, ^ 1 sh(y/2) + 3e^'' sh(3y/2) 

2 ch(y/2)+e2rf ch(3y/2) ' ^ ' 

As a check, one finds at low T two different 2-level type responses depending 
on the sign of d, namely, (m) ~ (3/2) th (3y/2) for positive D and (m) ~ 
(1/2) th (y/2) for negative, with the right spin values, 3/2 and 1/2 respectively. 

One more derivative gives the longitudinal fluctuations entering in the sus- 
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Table 2: 5 = 3/2 energy levels — /3e,„ — dm? +ym, ladder factors £^ = 15/4 — 
m(m+l) and transition frequencies PAmm+i = /3(em — £m+i) = (i(2m + 1) + ?;. 
Here 5(5 + 1) = 15/4 and ■^S{S + 1) = 5/4, and to spare exponentials we have 
shifted all energies by d{l/2)'^, that is /Se^ := + 1^/4- 



ceptibility X|| = P (m?) — {m)'^: 

/3 l + 2e^''(4chy+ ch2y) + 9e^'' rc_,/o^ /ois^ 

^""4 [ch(y/2)+c2'ich(3y/2)]2 ' I* - -^/^i • 

Consistently with the above magnetizations, d 3> 1 gives X|| — /?(3/2)^/ ch^(3?//2) 
while d ^ — 1 yields x\\ — Pi^/'^Y I ch^(y/2) as low-temperature 2-level asymp- 
totics. 

ii) transverse x_l(5 = 3/2). We next assemble the transverse suscepti- 
bility from the quantities of Table [2] 

XI. - :r^[3e-^/2i^(2d-y)+4e"^/2A'(+2/)+3e+^/2/^(2d+y)] , (5 = 3/2) . 

(3.19) 

with the custom K{X) = {c^ — \)/X. We have played again with the detailed- 
balance property K{—X) ~ cT-^ K{X) to arrive at a compact form [e.g., note 
that the middle term is either e-y/"^ K{+y) = e+^'Z^ K{-y) = sh {y/2)/{y/2)]. 

iii) temperature dependence of x\\ and x± for 5 = 3/2. In the 

unbiased 3^=0 case, Eqs. (|3.18p and p.l9p reduce to 

Both approach the Curie curve Xc — 5 /AT when d = D/T 0, as seen in 
the high-temperature range of Fig. [2] (now 5(5 + 1) = 15/4). Again, at low 
temperature x\\ has to catch up with the 2-levcl asymptote X|| — S'^/T (with 
5 = 3/2 for D > 0), increasing faster than l/T in the crossover range. 

The transverse x± tends to a constant value as T — > 0, but it displays a 
small and broad maximum around T ~ D/2, which was absent in 5 = 1 (and 
hardly visible here). The maximum can be obtained from the alternative form 
X± = [thd + I d (1 — the?)], where d{l — thd) adds a small bump to the 
monotonous thd as T decreases {d = D/T). 
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We close with the susceptibihty for the ensemble with anisotropy axes ori- 
ented at random. Figure [2] shows that X = ^X|| + fxi matches the isotropic 
Curie curve over a wider T range than X|| and x_l; actuaUy down to T ~ D- 
2D, where it eventuaUy deviates downwards. This behavior was also shown by 
S* = 1, and will be addressed in the next section for larger S values. 

4 superparamagnets 

Discussing the T dependence of the susceptibility of anisotropic spins S = \ and 
5* 3/2 we have come across some features of the curves that we would like to 
study more systematically, for several S. We want to check if those behaviors 
are specific of some spin values, how they evolve with 5*, and whether they 
survive the classical/continuum limit 5 — > oo. 

We have in mind three features: (1) the crossover of x\\ from the Curie law 
= ^I3S{S + 1) toward the 2-state response as T is lowered, (2) the peak in 
the transverse X±(.T) around the anisotropy constant D, and (3) the deviation 
of the orientationally randomized % from Curie at low temperature^ 

Technically, the study as a function of S is eased by the compact expressions 
of Sec. [21 We just need an algorithm building a table hke those of Sec. [Sj for a 
given S, and feeding the formulas for X|| and XJ- [Eqs. (|2.25|) and ()2.26p ] with 
the entries £„, i'^ = S{S + 1) - m{m + 1), and Amm+i = £m~ £m+i- 

4.1 scaling with S 

To compare results for different S meaningfully, we must specify which param- 
eters are kept constant, or how they are scaled with S. The same applies to the 
classical limit. This kind of specifications is needed in any sensible limit-taking 
procedure in physics (thermodynamical limits, continuum limit from mechan- 
ics, or from lattice discretizations, etc.). Different specifications/scalings give 
different results, possibly answering different questions. 

We will use the following prescription. We compare spins with different 

5 but having the same (maximum) energy. For the spin Hamiltonian TYo = 
—D — B ■ S, this entails to keep fixed 

scaling : D = const. B S ^ const. , (4-1) 

as the spin is varied. Then the energy differences between adjacent levels A ~ 
D S + Bz will decrease as A ~ 1/5, and the levels will approach each other 
accordingly. In Fig. [T] we have shown three spins with the same anisotropy 
barrier AU = DS^ but different number of levels (3, 5 and 9). Therefore, with 
this convention the question we will be answering is how the number of levels, as 
a measure of discreteness/quantumness, would affect the properties observed. 

In the rest of this section we will address the three points mentioned for the 
uniaxial model TLq — ~D — B^ Sz- Having in mind actual superparamagnets, 

One may think that (1) and (3) are the same, but in classical superparamagnets there is 
crossover in X|| but no deviation of ^ from Curie |39) . 
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we use D > and set = for simplicity. The quantity DS"^, being fixed, 
provides a natural scale of temperature T/DS^, with D ^ 10^-10'^ K in 
many superparamagnets O [TU] . Similarly, the output susceptibility curves will 
be normalized as x/S{S + 1), to facilitate comparisons. For example, all curves 
Tx would go over 1/3 at high T, and we can discuss how the low T behavior is 
modified by the discreteness of the spectrum. 




Figure 3: Longitudinal susceptibility x\\ vs. T for several S. Curves presented as 
Tx/S{S+l), to show the crossover from the high-T Curie regime, Tx/S{S+l) = 
1/3 (thin dashed line), to the 2-level regime at low temperatures. The arrow 
at ~ 1/25 indicates the lower limit of the observable equilibrium temperature 
window due to the onset of non-equilibrium effects (superparamagnetic blocking; 
cf. gOl). 



4.2 longitudinal response for various 5" 

As already mentioned, understanding well the T dependence of the susceptibility 
of the paramagnetic "ideal gas" is important because x(^) curves are routinely 
used as indicators of spin-spin interactions, anisotropies, etc P^l [5D1 15^ . 

In Sec. Owe saw that for anisotropic S* = 1 and S = 3/2, at high enough 
temperature, the susceptibility approaches the Curie asymptote Xc = ■^PS{S + 
1), whereas in the low T range the highly anisotropic 2-state regime xi = 
emerges (if you prefer, Heisenberg to Ising crossover). This is expected VS', 
because at high enough T the anisotropy term plays a minor role e~^^"^ ~ 1, 
and the spin effectively becomes a free (quantum) rotor. Indeed, arithmetic 
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sums like J2k=i = l^iN + 1){N + i) give 

s s 



J2 m'c-^'^ ^ ^\S{S + l){2S+l) , (4.2) 

rn— — S m— — S 

which divided by Zq 2± J2m 1 = 25* + 1 (number of states) produces the Curie 
law Xc ~ ■^f3S{S + 1), for aU S. Then, as ~(3em — (iDm^, the onset of deviation 
from isotropy can be defined as the temperature where e~^'^" = 1 and e~^^^ 
(population of the "poles" and the "equator") start to differ appreciably. Say 
by 5%, then e^^'/"^ 1 + DS^/T - 1.05 gives T/DS^ - 20. This estimate 
holds y S, as seen in the high-T range of Fig. [3] 

In the opposite low T regime, we can approximate X|| by using the lowest 
levels m = ±5 (effective 2-state; classically only the poles populated) 

(4.3) 

In this low T range, Tx\\ shows a plateau (Fig. [3]) which does not bend until 
the next levels e±(s_i) become appreciably populated. This suggests that the 
leaving this "Ising plateau" is governed by Tj := Ag^s-i = D{2S — 1) (the 
energy difference to the first excited pair). As this energy enters exponentially, 
the curves' bend would be negligible at a T five times lower [exp(— x) ~ at 
X ~ 4-6]. Indeed in our scaled units Ti/bDS^ = {2S - 1)/5S^ ~ 2/55. This 
gives a longer plateau the smaller S is, in agreement with the curves of Fig. [3l 
Note that for large S, not only is the Ising regime shorter, but it is reached 
more slowly when decreasing T (non exponentially)!^ 

With the two estimates above we have characterized the extent of the tem- 
perature crossover for all S. As for the magnitude of the "excursion" , it just 
follows from Tx\\ evolving from ^S{S + 1) up to 

5^ 35 f3 S ^ oo , , 

(4.4) 



i5(5+l) 5+1 13/2 5 = 1 • 

This is maximum classically, 3, and decreases with 5 to get halved for 5 = 1. 
This merely refiects the familiar quantum-mechanical fact of 5^ not having 
length 52, but 5(5 + 1). 

Having characterized the crossover, one could assess the temperature ranges 
where approximate modelizations can be employed. For example, for classi- 
cal nanoparticles the use of the Curie susceptibility is widespread (with the 
corresponding Langevin magnetization) , unfortunately well down to the super- 
paramagnetic blocking T / D ^ 0.04-0.1 (arrow in Fig. [3]). Here we see once 



For classical superparamagnets x\\ = P^i^R' /R where R{ot) = J (|dz exp{a2^) is the 
"partition function" with a = l^VjT, and R! = dR/do [9l l40l lill Ii2l I43lf44l . At low T (large 
ct) one can approximate R' / R ~ 1 — 1/a showing the non-exponential, power-law approach 
to X|| ^ At high T, on the other hand, letting a — > one has R' /R = /g dzz^ = 1/3^ 

recovering the classical Curie susceptibility xil — hP^i^ ■ 
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more that if a rough model is to be chosen, the 2-state model is preferable 
[with the associated (to) oc th (/iS/T)], as the pioneers of superparamagnetism 
properly did [1]. 

4.3 transverse susceptibility for various S 

The transverse response is a process quite different from the longitudinal one. 
X|| is about applying a small probe that shifts up and down the energy levels, 
and the associated re-population involving exp(— /3AC/) factors. However, a 
transverse probe mixes and splits the degenerate levels, and our intuition of the 
response gets somewhat lost. 

A classical cartoon can be of some assistance. In it the effect of b± is trying 
to rotate S out of the stable anisotropy minima (the "up" and "down" poles) 
toward some point in the equator. The result of this torque can be expressed as 
X± ^ \d^E. Thus, the transverse susceptibility is more sensitive to other fea- 
tures, like energy-well curvatures, rather than to barrier heights, which renders 
it as a valuable tool too. 

The torque picture holds at T = (Stoncr-Wohlfarth .2'_\). and we actually 
saw that, down there, X-L tends to some constant = D for S* = 1 and 
X_L = 3/41? for S = 3/2; Fig. d]). On the other hand, at high T > D the 
anisotropy potential becomes irrelevant and one should regain x± = ^I3S{S + 1) 
once again. If we insist in a Tx plot (Fig. [5]), this quantity will evolve from 
\S{S + 1) at high T down to 0, with an initial linear slope if Xj_(r) is nearly 
constant at low temperature. 

This behavior is what wc found before for S = 1 and S = 3/2. i.e., for 
the 3-lcvel and 4-lcvel systems. But 5 = 3/2 also exhibited a small and broad 
maximum in X-i-iT), which is absent in = 1. Is this a parity, integcr/half- 
integer effect? or the absence of peak is just an oddity of S" = 1 and it persists 
VS* > 1. The latter is suggested by the classical model having a maximum too. 
And indeed, plotting the bare x± vs. T (lower panel of Fig. 2]) we see that the 
peak is there for all > 1, becoming more apparent for large S. 

Some interpretation can be provided combining the torque picture above 
with thermal activation. A small T ^ can assist in leaving the potential 

The folk view associates superparamagnetism with energy barriers AU = D larger 
than T. But a crude estimate of the over-barrier relaxation time r = tq exp(/3AC/) gives a 
pre-factor tq ~ 10^'^-10~®s for molecular clusters and tq ~ 10~^"-10~^^ s for nanoparticles 
1101 15]. Superparamagnetism is to be observed for measurement times tm 2> t. Then one 
can have an equilibrium/superparamagnetic temperature range as wide as 25 > /3Af/ > 0, for 
static measurements tm ~ 1-100 s, showing that the folk ascription 1 > /BAU > is hopelessly 
restrictive. For example, almost the whole crossover from /3AU <S 1 to fSAU 2> 1 can fit in 
the observational equilibrium temperature window (we marked its lower limit ~ 1/25 in the 
T axis in Fig.^. 

The initial decrease, and then increase in x±(^ = 0), is due to our normalization and 
scalings (Fig. H bottom). For fixed D, the bare x±(S,T = 0) = 1/[2D(1 - 1/25)] [Eq. iji^S)! ] 
decreases monotonically with S: = ^/D, Xx(3/2) = 3/4D, Xx(2) = 2/SD, to get 

halved as x±{S oo) = 1/2D. But the scaling D = 1/5^, plus x±/S{S + 1) gives x± ^ 
S^/(2S^ -1-5—1), which starts from 1/2 at 5 = 1, decreases sharply to the minimum 4/9 at 
5 = 2, and slowly returns to 1/2 as 5 — > oo. Notice Xx(3/2) = x±(3) = 0.45 in these units. 
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Figure 4: Upper panel: Transverse susceptibility vs. T for various S, plotted as 
Txl S{S + 1). The two lowest curves correspond to those of S" = 1 and S = 3/2 
in Fig. [21 In the lower panel we plotted the raw susceptibilities, showing the 
peaks developed for S > 3/2, and how they evolve toward the classical curve. 

minima and, on average, may help reorienting S toward the transverse field 
(increasing the response). Too high a T, however, and the custom thermal 
misalignment would set in, decreasing the susceptibility. Then a peak in between 
seems natural from the competition of both processes. 

Well, but why is there no peak for 5 = 1? We may answer that, lacking 
intermediate levels, the T that assists rotation toward m = becomes too large. 
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and we only find the thermal decrease. Quantum mechanically it can be put in 
the following way. A transverse field bx mixes the states |m) and |m ± 1), and 
their contribution to (Sx) involves the factor \{m + l\Sx\'m)\'^ /{sm+i — £m) x 
(•g-/3em+i _ Q-I3^my Then the contribution of m = 5 — 1 is larger than that of 
m = 5, as the above factor increases when m climbs the ladder out of to = S* 
(mostly because the levels get closer if = —Drn^ so that — £„) 

becomes larger). However, for S" = 1. the level m = has no level above it 
providing such an increase. Then populating to = thermally does not increase 
X±^ and there is no peak for spin one. 

As for the behavior of the peaks with S", we see that they move only a 
little in the scaled units T/DS'^. This rules out a characteristic temperature 
of the type Tq £'(25' — 1) (as invoked in the longitudinal case), because then 
Tq/DS"^ ~ l/S, and the peak would shift left and disappear classically. We 
need another governing energy scale; the total barrier AU = D would do, 
since is constant in our units. Indeed, exp(— ACZ/To) ^ would mark the onset 
of misalignment, as the equatorial levels then start to be populated. Using again 
the rulc-of-thumb exp(— a;) ~ at a; ~ 4-6, we would have AU/Tq ~ 4-6. Then 
To/DS'^ - 0.16-0.25, compatible with the peaks' location in Fig.gl 

4.4 response for randomly oriented axes 

We conclude with the susceptibility of the ensemble with anisotropy axes dis- 
tributed at random. This case is of experimental interest in powdered samples 
and liquids (ferro-fluids [S]), as well as in solid systems with oricntational dis- 
orders. 

4.4.1 restoring the isotropy and the classical limit 

In classical superparamagncts the oricntational average 6| = 1/3 and b\~2/2> 
leads to a full restoration of isotropy [46l Fig. 1] [47] 

That is Ria) = /g^dz exp(az^), with a = AU/T, and its a-dcrivative R' dis- 
appear from X- This washing away of any trace of the anisotropy constants is 
not specific of the uniaxial model TCq ~ ~D 5^, but it holds for any classical Ho 
with the reflection symmetry Ho{—S) = Ti-oiS) (see [31] and references therein). 

Quantum mechanically, however, we know that this cannot be exact. In 
Fig. [2] we showed that for S — 1 and S = 3/2 the random axes averaging indeed 
restores isotropy, x — Xc7 over a wider temperature range. And that this occurs 

In the classical model the peak can be found plotting x± = PtJ-^iR — R')/'iR [9], and is 
located around To ~ 0.23 1401 . On the other hand, the effective eigenvalue method 1451 gives 
an effective relaxation time t^j showing a similar maximum [9l p. 126]. The reason seems to 
be that Tef can be expressed in terms of equilibrium averages, and happens to be proportional 
to 1 - (z^) oc XX. 
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Figure 5: Random axes susceptibility vs. temperature for various S, presented 
as Tx/ S{S + 1). (Notice the factor 1/2 difference for 5 = 1 with respect to 
Fig. [21 due to the normahzation of The fowcr panel shows the same curves 
but over a wider temperature range in logarithmic scale, to test the approximate 
formula (|4.7|) (small circles). 



even with X|| well on its way to xi = PS'^j but at a certain temperature the 
susceptibility starts to deviate noticeably downward (it cannot catch up with 
^(38(3 +1)). That temperature is of the order of £> for S* = 1 and 2D for 
S ~ 3/2 (Fig. (2); suggesting that the relevant/governing energy scale is again 
Ti ~ D{2S — 1). Indeed, in units of D the extent of the bent range would be 
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Table 3: Lowest energy levels for spin S at zero field —(3 Em = dm?, lad- 
der factors £^ = S{S + 1) — m{m + 1), and the transition frequencies 
PAmm+i = /3(£m " £m+i) = d{2m + 1). The approximate partition function is 
~2(e-'3^s +e-'3^s-i) =2e-'^'=s(i + e-") (cf. gi App. A]). 

from T = to Ti/D S"^ ^2/S (see Fig. [S]). This smoothly gives the connection 
with the classical result of full restoration of isotropy by letting \/ S ^ 0. 

4.4.2 approximate formula for % 

To confirm the previous estimates and scalings we have derived an approximate 
formula for x based on a few-levels treatment. It may also provide some insight 
in the origin of the deviations from Xc = \PS{S -{- 1). 

We started looking at Eq. (|2.26p for xj.- One notices that letting K„i 
1 (for instance, considering K{PA) plus A ~ l/S), the susceptibility is left 
as xr = {P/2){el), with (£1) := Em^le-^^'- /Zo. Then, combining the 
average of the ladder factor £^ = S{S -t- 1) — m{m + 1) with the longitudinal 
susceptibility gives |x|| +|x!l° = ^0S{S + 1), explaining the notation x±°- The 
actual susceptibility is obtained by adding and subtracting (i'^) to Eq. (|2.26p 

X±^xT + iP/2){iliK^-l)) . (4.6) 

This form, times 2/3 and combined with the longitudinal part (|2.25p . gives 
X = Xc + (/3/3) {i-f„ {K„i — 1)). Thus in the second term (of transverse origin) 
we have isolated the source of deviation from isotropy, but we have not done 
approximations vetF^ 

To compute the suspect term (^i'^{Km — 1)) we just include the contribu- 
tion of the two lowest levels (at each side; see Table [3]), getting the following 
susceptibility 

^4e+ Li.(l-l/25)!;e-n ' ^-P^S-l). (4.7) 
^ ^ ^ 

By analogy with the classical superparamagnetic susceptibilities l|4.5|l one could introduce 
{rn?) =: S{S + 1) R' /R (which defines R' /R), and write (f = T/S{S + 1)) 

fx\\ =R'/R f x± = (■R--R')/2if. + corrections 

The first term in x± corresponds to Xx° (i-C-i when combined with X|| and randomized gives 
Curie), and the corrections correspond to (^(.'^{Km — !))• 
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From this expression we identify the Ising x\\ — PS"^ at low temperatures. Be- 
sides, invoking U := (3D{2S — 1) ^ 1 at low T, one can read off the torque 
transverse susceptibility V5, namely 

This gives the familiar values x± = 1/-D for S* = 1 and x± = 3 /AD for 5* = 3/2. 

Equation (|4.7p is governed by the level difference fl := /3A5 5_i, as expected 
from our lowest-level approximation. We have plotted it together with the exact 
X curves in the lower panel of Fig. [5l The description it provides is reasonably 
good at all temperatures, for all spin values, and quite good for 5* = 3/2 and 2. 

Multiplying now Eq. ((47)) across by T and using - 1/25") = n/2S we 

can write the compact form 

Sf 2 1-C- " 
3 I ^J^TTi" 



Tx^-[S+--rT^] , n = PD{2S-l). (4.9) 



This gives the S'-dcpcndcnt intercept of the x a-xis at T = of Fig. O and the 
initial linear growth at low temperatures (1/ri oc T with c^^^ ~ 0). But Eq. (|4.9p 
happens to capture as well the restoring of Curie at high T (fi <C 1). Indeed, 
Taylor expanding (2/17)(l - e"")/(l + e~") ~ 1 - 17^12, the constant term 
nicely produces ^S{S + 1), the Curie constant. The second term can then be 
used to estimate the onset of deviations from Xc- Fixing a 5% deviation, i.e., 
~ 0.05, and writing Q = {DS'^/T){2S - l)/S^, yields T/DS^ - 5/25, 
in agreement with our previous estimate 2/S. 

As for the use of the approximate x one should bear in mind that Eq. (|4.7p 
was derived under low T conditions (few populated levels). However, the reason- 
able agreement with the exact curves in the whole temperature range, and for 
all S, indicates that it could be used safely in modelization of uniaxial magnets 
with axes distributed at random. 



5 summary 

The understanding of the properties of paramagnets belongs to a long tradi- 
tion linking magnetism, quantum mechanics and statistical mechanics. Our 
aim here has been to extend the theoretical framework to permit the study of 
some equilibrium problems for arbitrary values of the spin. In this frame, one 
can connect from landmark results for quantum paramagnets (Curie-Brillouin, 
transverse response of anisotropic spins) all the way up to the theory of classical 
superparamagnets, developed for magnetic nanoparticles, and revived with the 
young molecular magnetic clusters (single-molecule magnets) . We have used the 
language of magnetism throughout, but the formalism is closely related with the 
effective big-spin description of collections of 2-level systems used in atom optics 
and two-mode Bose condensates [191 HQI [21] . 

We focused on uniaxial spins and the temperature dependence of the mag- 
netic susceptibility, due to its traditional significance and its routine use as 
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characterization tool. We investigated three features: (1) the crossover of the 
longitudinal susceptibility from Curie to the 2-state regime, induced by the mag- 
netic anisotropy, (2) the peak in the transverse x± vs. T, and (3) deviations of 
the random ensemble susceptibility from the Curie law (absent in the classical 
limit). We identified and characterized the relevant energy/temperature param- 
eters governing the phenomenology, and studied how they scale with S. We did 
this starting from small spins 5 = 1, 3/2, 2, then moderate S' = 5, 10, . . . and 
eventually big spins S = 30, 50, 100, connecting with classical superparamag- 
netic phenomenology. 

The equations we employed do not rely on Van Vleck's method, as they 
follow from the general Kubo correlator formalism of linear-response theory. 
We worked and particularized this formalism, with the above problems in mind, 
to produce ready-to-use formulas [Eqs. (|2.25[) and (|2.26p ]. which only require 
the input of the unperturbed spectrum and angular-momentum ladder factors. 

We also derived approximate expressions and assessed their ranges of valid- 
ity. This turned out to be quite good for the lowest-levels approximation (|4.7p 
to the susceptibility of the randomly-oriented ensemble. It could be safely used 
as a compact modelization of susceptibility data for arbitrary S in such an ex- 
perimentally relevant case. In real magnets, however, the unperturbed Hamil- 
tonian will include terms non-diagonal in the standard basis, like D S^. or KS^.. 
These terms were not accounted for in the derivation of the approximate ex- 
pression (|4.7p and, when relevant, one should turn back to the more general 
equation ((2^ . 
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